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Within the framework of one-dimensional Laval-Dubrulle-Nazarenko type model for the La- 
grangian acceleration in developed turbulence studied in the work [A.K. Aringazin and M.I. Mazhi- 
tov, cond-mat/0305186 we focus on the effect of correlation between the multiplicative noise and 
the additive one which models the relationship between the stretching and vorticity, and can be seen 
as a skewness of the probability density function of some acceleration component. The skewness of 
the acceleration distribution in the laboratory frame of reference should be zero in the ideal case 
of statistically homogeneous and isotropic developed turbulent flows but when considering acceler- 
ation component aligned to fluid particle trajectories it is of much importance in understanding of 
the cascade picture in the three-dimensional turbulence related to Kolmogorov four-fifths law. We 
illustrate the effect of nonzero cross correlation parameter A. With A = —0.005 the transverse (x) 
acceleration probability density function turns out to be in good agreement with the recent exper- 
imental data by Mordant, Crawford, and Bodenschatz. In the Random Intensity of Noise (RIN) 
approach, we study the conditional probability density function and conditional mean acceleration 
assuming the additive noise intensity to be dependent on velocity fluctutions. 

PACS numbers: 05.20.Jj, 47.27.Jv 



I. INTRODUCTION 

Tsallis nonextensive statistics [l| inspired approach |3, 
0, was recently used 0, @ to describe Lagrangian ac- 
celeration of fluid particle in developed turbulence; see 
also 0) B 13 • I n Ref. ITcjl w e reviewed some refinements 
of this approach 0, Il2lll3| . 

Review and critical analysis of the applications of var- 
ious recent nonextensive statistics based models to the 
turbulence have been made by Gotoh and Kraichnan . 
An emphasis was made that some models lack justifica- 
tion of a fit from turbulence dynamics although being 
able to reproduce experimental data to more or less ac- 
curacy. A deductive support from the three-dimensional 
Navier- Stokes equation was stressed to be essential for 
the fitting procedure to be considered meaningful. 

Recently Laval, Dubrulle, and Nazarcnko [lj| have de- 
veloped a stochastic kind of Batchelor-Proudman rapid 
distortion theory approach to the three-dimensional 
Navier-Stokes equation using separation of large-scale 
and small-scale velocities and Gabor transformation (lo- 
calized wave-packets) to derive one-dimensional Langevin 
toy model for small-scale velocity increments both in the 
equally interesting Eulerian and Lagrangian frames; see 
also recent paper 0. The large-scale terms entering 
the resulting approximate small-scale equation are re- 
lated to large-scale strain and inter-scale coupling and 
are treated as noises with a given statistics. The small 
scales are stochastically distorted in certain way as a com- 
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bined effect of the large scales and the inter-scale cou- 
pling. Short-time correlated character of the distortions 
follows from the numerical study of decaying turbulence; 
longtime correlations and dependence on velocities which 
may be present here as well have not been modeled in 
the first step. This approach allows one to account for 
nonlocal interaction effects in small-scale turbulence via 
a simple random multiplicative process driven by cou- 
pled Gaussian white-in-time multiplicative and additive 
noises, while local small-scale interactions are modeled 
by a turbulent viscosity. 

In a comparative analysis of some recent one- 
dimensional Langevin toy models of fluid particle accel- 
eration in developed three-dimensional turbulence 01 
we have demonstrated that the one-dimensional Laval- 
Dubrulle-Nazarenko (LDN) type model 0, 0] , with the 
model turbulent viscosity v t and delta-correlated Gaus- 
sian white multiplicative and additive noises, formulated 
for the Lagrangian acceleration meets the experimental 
data on acceleration statistics ^3,^3 to a good accuracy. 
Particularly, it was shown that the resulting contribution 
to fourth order moment, a 4 P(a), does peak at the same 
values as the experimental curve, in contrast to predic- 
tions of the most of other stochastic models [jj El M ^2 • 

Also, within the framework of Random Intensity of 
Noise (RIN) approach 0| to the LDN type model the 
assumption that the additive noise intensity a depends 
on absolute value of velocity fluctuations u was found to 
imply the conditional probability density function P(a\u) 
and the conditional acceleration variance {a 2 \u) which 
are in a good qualitative agreement with the recent ex- 
perimental data on the conditional acceleration statistics 
reported by Mordant, Crawford, and Bodenschatz [l9| . 
These results have been obtained in the particular case 
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when the correlation between the multiplicative and ad- 
ditive noises is taken to be zero. 

In the present paper, we fill the gap by studying the 
effect of nonzero cross correlation of the noises that mod- 
els a relationship between stretching and vorticity in the 
three-dimensional case ^| and can be seen as a skew- 
ness of the probability density function of acceleration 
component. 

We remind that shell models of turbulence arc in- 
capable to describe the observed skewness generation 
along the scale of the probability density function of 
the longitudinal Eulerian velocity increments [l5j|. This 
skewness is of much importance in understanding of the 
Richardson-Kolmogorov cascade picture of the homoge- 
neous isotropic turbulence. The skewness is particularly 
related to a non-zero value of the Eulerian third-order 
velocity structure function, and therefore to the essence 
of the turbulent cascade via the Kolmogorov four-fifths 

law [THHEI. 

Intermittency of homogeneous and isotropic turbu- 
lence is usually characterized by a nonlinear dependence 
of the scaling exponents on the order n of moment, and is 
well established experimentally both in the Eulerian and 
Lagrangian frameworks for approximately homogeneous 
and isotropic flows. 

The Eulerian and Lagrangian anomalous scalings trace 
back to a local inhomogeneity of the flow and long- 
time correlations in the particle accelerations respec- 
tively [22|. The probability density functions of the Eu- 
lerian and Lagrangian velocity increments both exhibit 
the same behavior: they are approximately Gaussian at 
large scales and progressively develop stretched exponen- 
tial type tails when the spatial and time increments de- 
crease down to the Kolmogorov length and time, respec- 
tively. In the limit of zero time increments the proba- 
bility density function of the Lagrangian velocity incre- 
ments converges to that of the Lagrangian acceleration. 
In practice, one uses nonzero time increments r lying 
within the range of strong viscous dissipation, such that 
velocities are smoothed and the characteristic relation 
u(t + t) — u{t) = ra(t) holds to a good accuracy. This 
time scale is known to be about or less than the Kol- 
mogorov time. 

In the Eulerian framework, the relative scaling expo- 
nents of the absolute moments of the longitudinal Eu- 
lerian velocity increments, (\u(x + r) — u(x)\ n ), were 
measured to be Of = 0.36, 0.70, 1.28, 1.53 for n = 
1, 2, 4, 5 [2^, and the third-order Eulerian velocity struc- 
ture function, taken as a reference, is known to scale lin- 
early with the spatial separation r due to the Kolmogorov 
scaling; see, however, Ref. |2l|. in which Kolmogorov- 
Novikov equation accounting for the external forcing 
length in the r 2 term has been confirmed experimentally 
for an impressive range of Reynolds numbers to a very 
high accuracy. 

This picture converts into the Lagrangian framework. 
The relative scaling exponents of the absolute moments 
of the Lagrangian velocity increments (one component), 



(\u(t + r)-u(t)\ n ), were measured to be C„ = 0.56±0.01, 
1.34±0.02, 1.56±0.06, 1.8±0.2 for n = 1,3,4,5 ||, and 
the second-order Lagrangian velocity structure function, 
taken as a reference, is known to scale linearly with 
the time increment r due to the dimensional analysis. 
One therefore expects skewness of the probability density 
function of the Lagrangian velocity increments in time 
since ^ is nonzero, and hence that of the Lagrangian 
acceleration. 

The one-dimensional LDN toy model was formulated 
originally in both the Eulerian and Lagrangian frame- 
works for velocity increments in the frame comoving with 
the wavepacket [l^. We use the Lagrangian formula- 
tion and the exact result for probability density func- 
tion for LDN type model obtained as a stationary so- 
lution of the Fokker-Planck equation associated to the 
one-dimensional Langevin equation for the acceleration 
component [Iclll5j . 

d t a= {£-v t k 2 )a + a ± . (1) 

Here, the Gaussian white noises £ and <r±_ model stochas- 
tic forces in the Lagrangian frame and are defined by 

<£(*)> = 0, mt(t'))=2D5(t-t'), 
(a±(jt)) = 0, (a±(jt)ax(^)) = 2aS(t - t'), (2) 
(m<r±(t')) = 2X6(t-t'). 

All the noises are treated along a particle trajectory. 
Here, the free parameters a, D, and A measure intensities 
of the noises and their cross correlation, respectively. 

Such a choice of noises is motivated not only by sim- 
plicity of their statistics but also by the DNS of decaying 
turbulence |l5j . 

Similar problem within the Eulerian framework has 
been recently investigated, with the result that the noise 
entering Langevin type equation can be safely taken 
delta-correlated in scale. Namely, the Eulerian exper- 
imental study of Langevin modeling of velocity incre- 
ments by Renner, Friedrich and Peinke |24| and Marcq 
and Naert [2(j reveals that the longitudinal velocity is 
correlated over distances much larger than the correla- 
tion length of its spatial derivative, so that the Markovian 
approximation is accurate in the inertial range. Approxi- 
mation of a short-correlated noise by the delta-correlated 
one is usually made due to the scale hierarchy validating 
the use of Langevin type equations. It is natural to map 
this result to the Lagrangian domain. 

In the Eulerian LDN framework, the cross correlation 
parameter was found to control third-order longitudinal 
velocity structure function due to a kind of generalized 
Karman-Hovarth relationship derived in Ref. 15]. Since 
the noise distributions are taken to be not skewed, this 
suggests that the parameter A defined in Eq. J2J should 
be nonzero. 

Experiments |2(j show that the distribution of the Eu- 
lerian longitudinal velocity increments is slightly asym- 
metric in the inertial range of scales, with the skewness 
factor being about S — —0.25 for the studied R\ = 430 
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flow; the odd-order moments (2n + 1) are known to be 
small as compared with even-order ones (2n + 2). The 
Langevin model of Ref. [2£| assumes the use of only one 
noise driving Eulerian velocity increments across scales, 
and this noise is characterized by a slightly skewed distri- 
bution (the skewness factor S — 0.55) and long tails (the 
flatness factor F = 8.5 as compared to F — 3 for a Gaus- 
sian). Gaussian approximation for this noise was used 
to derive the Fokker-Planck equation. It was found that 
this approximation and accounting for corrections com- 
ing from non-Gaussianity of the noise imply some persis- 
tent deviations from the observed scaling of third-order 
and fifth-order velocity structure functions. This may 
indicate that some different type of Langevin or Fokker- 
Planck equation should be used as ansatz. 

The probability density function as a stationary solu- 
tion of the Fokker-Planck equation associated to Eqs. JTJ 
and J2J was calculated exactly 0] , 

CeM-Vtk 2 /D + F(c) + F(-c)} 
{a> {Da 2 -2Xa+a) 1 / 2 {2Bka + Vt k 2 ) 2BXk / D2 ' 1 ' 



where we have denoted 



F(c) = 



ci 



k 2 



ln[- 



2D : ' 



2c 2 D 2 c cic 2 (c-Da + \) 
x (B 2 (X 2 + cX- Da)a + c{Dv 2 t k 2 + c 2 v t ))], (4) 

c = -i^Da-X 2 , v t = ^Jv 2 + B 2 a 2 /k 2 , (5) 
ci = B 2 (AX 3 +AcX 2 -3DaX-cDa) + D 2 (c+X)^k 2 , (6) 



c 2 = JB 2 {2X 2 + 2cX - Da)k 2 + D 2 v 2 k A , (7) 



and C is normalization constant. 

Without loss of generality one can put k 



1 and a — 1 



by appropriate rescaling of the parameters D, B, Vq, and 
A yi| to make a fit of the above P(a) to the experimental 
data. As one can see from Eq. 10, the parameter A 
introduced in Eq. (0) is responsible for an asymmetry of 
the distribution with respect to a — > —a. 

The fit for the particular (symmetric) case, A = 0, has 
been made in Ref. 01 ■ The flatness factor of the stud- 
ied distribution which characterizes the widening of its 
tails (when compared with a Gaussian) is found to be 
F = 42.5 (for fc = 1, a = 1, D = 1.130, B = 0.163, 

= 2.631, C = 1.805) that deviates from the flatness of 
the experimental curve, In numerical cal- 

culations we used a cutoff by restricting the integration 
range by \a\/{a 2 y^ 2 < 1000. 

The RIN approach extends the LDN type model 
by assuming certain relationship of noise intensities and 
in general other model parameters to velocity fluctua- 
tions u, and enables one to study acceleration statis- 
tics conditional on velocity fluctuations [Io|. A depen- 
dence of acceleration distribution on velocity fluctuations 
is known to violate Kolmogorov 1941 local homogeneity 
of the flow [3. 

The paper is organized as follows. 
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FIG. 1: Acceleration probability density function P(a). Dots: 
the experimental data for the transverse component of ac- 
celeration at R\ — 690 by Crawford, Mordant, and Bo- 
denschatz |T3| . Dashed line: stretched exponential fit |i~8j|. 
13 = 0.513, cr = 0.563, 7 = 1.600, C = 0.733. Solid line: The 
model k = 1, a = 1, D = 1.100, B = 0.155, u = 2.910, 
A = -0.005, G = 3.230. x = a/{a 2 ) 1/2 denotes normalized 
acceleration. 
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FIG. 2: The contribution to fourth order moment, a 4 P(a) 
Notation is the same as in Fig.0 



In Sec. [n] we study a skewness effect implied by a 
nonzero cross correlation parameter A, both for uncon- 
ditional and conditional probability density functions of 
the acceleration component. In Sec. II I II we study con- 
ditional mean acceleration using the RIN extension of 
the LDN type model ©. In Sec. llVl we summarize the 
obtained results and make a few concluding remarks. 



II. THE SKEWNESS 

A. The unconditional probability density function 

In the present paper we generalize the consideration 
made in Ref. |10| by letting the cross correlation param- 
eter to be nonzero, A 7^ 0. For this case, the exact prob- 
ability density function @ will be used. 
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A sample fit of the distribution J2J is presented in Fig.Q] 
and the contribution to the fourth order moment is plot- 
ted in Fig. [21 in which the signature of a skewness of the 
experimental distribution (dots) of the acceleration com- 
ponent can be seen as a small difference of about 5% in 
heights of the two peaks. Geometrically, this component 
of acceleration corresponds to the direction transverse to 
the axial symmetry (z) axis of the large-scale forcing of 
studied flow [l7| . 

We remind that the turbulence was generated in a flow 
between counter-rotating disks in a cylindrical container, 
and the flow significantly deviates from the ideal of homo- 
geneity and isotropy. By the large-scale flow symmetry 
the two transverse components (x and y) are taken sta- 
tistically equivalent, and distinct from the axial compo- 
nent (z), in the Cartesian laboratory frame of reference. 
Particle accelerations were measured in a small volume 
in the center of the flow chamber within which the flow 
is approximately homogeneous and isotropic. Only one 
transverse (x) and the axial (z) components were actu- 
ally measured. Statistical properties of the unmeasured 
y component are expected to be identical to those of the 
x component p^lllSj. Throughout the paper we consider 
the data for the x component of acceleration. 

One observes a better agreement of the skewed curve 
(solid line in Fig. |2J for the fitted value 

A = -0.005 (8) 

{k = 1, a = 1, D = 1.100, B = 0.155, v Q = 2.910, C = 
3.230) with the data points in the intermediate range 
of positive and negative accelerations, ^/(a 2 ) 1 / 2 ~ 10, 
as compared with (symmetric) curves implied by the 
stretched exponential fit (dashed line in Fig. |2J [l3 , the 
RIN chi-square Gaussian model fit the RIN log- 

normal model 0, and the Reynolds model 

We note that better result was recently obtainedJl8| 
by processing the original (unfitted) curve of Ref. ({j in 
exactly the same way as that for the data but some de- 
parture from the experimental data still persists. Par- 
ticularly, the associated contribution to fourth order mo- 
ment, a 4 P(a), does not peak at the same value as that 
for the experimental curve, being however very close to 
it. 

Since in statistically homogeneous and isotropic tur- 
bulent flows for the x, y, or z component of accelera- 
tion one should have zero skewness, we could attribute 
the observed small skewness to anisotropy of the studied 
Rx = 690 flow. 

We remind that the relative root-mean-square (rms) 
uncertainty of the experimental P(a) is about 3% for 
^/(a 2 ) 1 / 2 < 10 (the most accurate part of the distri- 
bution), and is less than 40% for \a\/(a 2 ) 1/2 < 40 [H|. 
Clearly, less uncertainties are required to obtain good 
quantitative description of the skewness as it appears to 
be a quite tiny effect. Nevertheless, in the present section 
we use the available experimental data to demonstrate 
the effect of nonzero A. 



Recently reported high precision data [19j show that 
the mean acceleration conditional on velocity fluctua- 
tions is nonzero and increases with the increase of ve- 
locity. This was claimed to be related to the anisotropy 
of the studied flow, although it was pointed out that DNS 
of homogeneous isotropic turbulence also shows slightly 
nonzero mean acceleration. We will consider this issue 
below and in Sec. IIIII 

It should be emphasized that for the studied R\ = 690 
flow the ratio between the variances of the x and z com- 
ponents of acceleration was measured to be slightly differ- 
ent from unity due to ao x /ao z ~ 1.06 |17|. It was pointed 
out that experimental biases do not affect acceleration 
measurements (they depend mostly on the velocity) ex- 
cept to the extent that the acceleration and velocity are 
correlated. The shapes of the probability distributions 
for the x and z components of acceleration were mea- 
sured to be approximately the same at high Reynolds 
numbers. This level of x-to-z anisotropy was found to 
persist for higher Reynolds numbers (data presented for 
R\ = 970) that may indicate a fundamental character of 
this phenomenon. Namely, while the K41 theory postu- 
lates complete universality it is still an open question to 
what extent statistical properties of the 3D turbulence in 
the inertial range do not depend on the details of large- 
scale forcing. As the presence of finite injection scale is 
felt through the entire inertial range via the anomalies 
of the Eulerian and Lagrangian scaling exponents (K62) 
it is natural to expect that the statistics of fluid parti- 
cle acceleration, which is generally associated to small 
scales of the flow, should reflect the anisotropy of the 
large scale forcing. From this point of view, the obser- 
vation supports the view that the induced anisotropy in 
the acceleration statistics is a rule rather than exception 
in the context of developed turbulence dynamics: Forced 
anisotropy at large scales is not washed out in the inertial 
range of a high-Reynolds-number flow and seems to be 
felt at small scales, smaller than the Kolmogorov length. 
This kind of anisotropy in acceleration statistics may not 
be directly related to the parameter A. 

The cross correlation parameter A is nonzero for lon- 
gitudinal and zero for transverse velocity increments by 
construction [l^. For the Lagrangian acceleration, this 
reads that for the component of Lagrangian acceleration, 
a r , pointed along the corresponding Lagrangian velocity 
at some point of the particle trajectory the parameter 
A is nonzero while for the transverse component a n it is 
zero. The experimental data on a x and a z time series do 
not allow one to extract separately a T and a n , in order to 
verify that the a T acceleration PDF is skewed while the 
a n acceleration PDF is symmetric relative to the change 
of sign of acceleration. This requires obtaining the data 
on the whole set of the acceleration and corresponding 
velocity components, to have an access to geometry of 
individual trajectories of the tracer particle. Also, one 
should suppress the effect of above mentioned induced 
anisotropy. The influence of the large-scale anisotropy 
on a T and a n acceleration PDFs may occur to be differ- 
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ent. In contrast to the experiment where it seems to be 
difficult to provide sufficiently high level of isotropy and 
extract statistical data on a T and a n , the DNS could be 
used to test skewness of the a T and a n acceleration PDFs. 
The isotropy is well satisfied in DNS and at least it is free 
from strong large-scale anisotropy of the experimentally 
studied von Karman flow. 

The obtained small value J8J of the fitted cross cor- 
relation parameter A, as compared with the used values 
of the noise intensities D and a, is in a good agreement 
with the results of numerical Rapid Distortion Theory 
(RDT) analysis of the noise cross correlators in decaying 
turbulence 15]. Particularly A turned out to be about 
two orders of magnitude smaller than D and a. The time 
scale of the cross correlation as well as of the autocorre- 
lation of additive noise. 

Below we use the RIN extension ^(j of the LDN model 
to obtain and study the conditional probability den- 
sity function. 



B. The conditional probability density function 

In the RIN approach the result is treated in gen- 
eral as a probability density function conditional on the 
parameters involved in the model. It was found |lfj| that 
for A = only a variation of the additive noise intensity 
a, with a = e u / uo , qualitatively meets (i) the observed 
variation of the shape of experimental probability den- 
sity function of the transverse component of acceleration 
conditional on the transverse component of velocity fluc- 
tuations u with variation of u, and (ii) the polynomial 
type increase of the normalized conditional acceleration 
variance with an increase of |u| |l9|. 

In the present paper we study the effect of nonzero A. 
In general one may expect that the variation of A causes 
not only essential variation of the skewness but at the 
same time a considerable variation of the acceleration 
variance. 

The experimental conditional distributions P(a\u) for 
u ranging from to 3 rms velocity were found to be 
almost of the same shape as the experimental uncondi- 
tional distribution P (a) 19]. This means that they share 
the same process underlying Lagrangian turbulence inter- 
mittency. 

In a qualitative analysis, we can use the fact that the 
experimental unconditional distribution P(a) is approxi- 
mately of the same form as the experimental conditional 
distribution P(a\u) at u — 0. As we suppose A = X(u) 
the net effect of a variation of u will be some variation 
of the parameter A in P(a). For this case, the result 
of our check is that the variation of the shape of curve 
with variation of A around the value A = —0.005 (the 
other parameters fixed) is not in a qualitative agreement 
with the experimental data on P(a\u) shown in Fig. 6a 
of Ref. Namely, the obtained plots shown in Figs. 01 
and ^demonstrate that the variation of A affects mainly 
the skewness rather than the variance of the distribution 
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FIG. 3: The acceleration probability density function. Dots: 
the experimental data for the transverse component of ac- 
celeration at i?^ = 690 by Crawford, Mordant, and Bo- 
denschatz |T3l . Dashed line: stretched exponential fit |l8j|. 
f3 = 0.513, a = 0.563, 7 = 1.600, G = 0.733. Solid lines: 
the model © at A = -0.005, -0.025, -0.05 (k = 1, a = 1, 
D = 1.100, B = 0.155, u = 2.910). x = a/(a 2 ) 1/2 . 
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FIG. 4: The contribution to fourth order moment, a 4 P(a). 
Notation is the same as in Fig. 01 

©■ 

We conclude that within the framework of RIN ap- 
proach the cross correlation parameter A (as well as the 
parameters D, i>q, and B as shown in Ref. 10]), with 
A = A(u), could not be responsible for the specific visible 
change of the shape of experimental conditional accel- 
eration probability density function with an increase of 
velocity fluctuations [Tsj . 

We are thus left with the only possibility: to assign 
the observed essential dependence of P{a\u) on velocity 
fluctuations u to the additive noise intensity a. In gen- 
eral, this is in an agreement with the RDT approach by 
Laval, Dubrulle, and Nazarenko, particularly with the 
approximate LDN model ||3J, in which only the additive 
noise intensity is characterized by a dependence on small 
scale velocity fluctuations coupled to large scale velocity 
fluctuations. 

It should be emphasized here that in the original LDN 
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the conditional mean acceleration only slightly deviates 
from zero, 



0.5 1 1.5 2 2.5 



FIG. 5: The conditional mean transverse acceleration as a 
function of standardized velocity fluctuations u for a = e' "' //3 ; 
k = 1, D = 1.100, B = 0.155, v = 2.910, A = -0.005. Tri- 
angles: normalized to unit conditional acceleration variance, 
(a\u) / (a 2 \u) ll/2 . Boxes: normalized to unit conditional accel- 
eration variance at u = 0, (a\u) / '{a 2 \Q) 1 ^ 2 . 



model the other parameters of the model, B, vq, and A, 
were not assumed to depend on small scale velocity fluc- 
tuations, and the multiplicative noise intensity D was 
found to depend only on large scale velocity fluctuations. 
Nevertheless, in Ref. [l0( we have evaluated the effect 
caused by each of these parameters, D(u), B(u), vo(u), 
and A(it), in order to verify independently whether this 
may qualitatively correspond to the set of experimen- 
tal conditional distributions P(a\u), u/{u 2 y/ 2 s [0,3]. 
Despite we obtained a negative result, presence of some 
weak dependencies of these parameters on u can not be 
ruled out on the basis of made qualitative comparison. 

Adopting the point of view that a — a(u) as a first 
approximation, in the next section we turn to a consid- 
eration of the mean acceleration conditional on velocity 
fluctuations u. 



III. THE CONDITIONAL MEAN 
ACCELERATION 

The conditional mean acceleration should be zero in 
homogeneous and isotropic turbulence, and departures 
from zero reflect the anisotropy of the studied flow al- 
though DNS of homogeneous isotropic turbulence has 
also shown slight departures from zero |19| . 

As the first step, the conditional mean acceleration can 
be calculated under the assumption that only a depends 
on velocity fluctuations u. We take an exponential de- 
pendence, a = e u / w °, uq = 3, which was found to be rel- 
evant from both the theoretical and experimental points 
of view [To| . 

Assuming that the unconditional distribution P{a) is 
approximately of the same shape as the conditional dis- 
tribution P(a\u) at u = 0, and using the same set of fitted 
parameters as in Sec. Ill Al (A = —0.005) we obtain that 



{a 2 \u) 1 / 2 



0.002, 



(9) 



for all values of |u| ranging from zero to 2.5 rms veloc- 
ity, with the tendency to decrease down to zero with the 
increase of as shown in Fig. (triangles). 

In general, this is in agreement with almost symmetric 
shapes of the experimental probability density functions 
of the component of acceleration conditional on the same 
component of velocity fluctuations [l^. Alas, an illus- 
trative character of the presented experimental plots and 
the increasing experimental uncertainty of P(a\u) at big 
■[^/(it 2 ) 1 / 2 does not allow us to make a definite conclu- 
sion since the skewness effect is very small to be readily 
seen from the experimental P(a\u). Also, we note that 
the predicted value (a|0)/(a 2 |0) 1/2 ~ 0.0024 (see Fis " 
does not contradict to that shown in Fig. 6b of Ref. [li 
(a|0)/(a 2 ) 1//2 < 0.05 (uncertainty of the presented data 
points does not allow us to give the upper bound more 
precisely) . 

The above result on the conditional mean transverse 
acceleration (Fig. [5j is however in a sharp contrast with 
the reported experimental dependence of the conditional 
mean transverse acceleration on it (Fig. 6b of Ref. 0) 
which displays that (a\u) / (a 2 ) 1 / 2 increases from about 
zero at \u\ — to about 0.3 at \u\/{u 2 y^ 2 = 2.5, in some 
nonlinear way. 

In essence the experimental data show that with the 
increase of velocity fluctuations |wj the mean of the con- 
ditional transverse acceleration distribution P(a\u) nor- 
malized to unit unconditional acceleration variance in- 
creases by small but appreciable amount. One would 
like to know whether this holds for the z component of 
acceleration. We expect that the mean for this compo- 
nent is bigger than that for the x component, partially 
due to about 3% smaller variance of the z component of 
acceleration |l7j . 

More detailed study is required to explain this highly 
remarkable phenomenon, which may indicate stronger 
coupling of the multiplicative noise to the additive one for 
bigger velocity fluctuations \u\ in the Lagrangian frame. 

Below, we use definitions of the noises provided by 
the LDN approach [l^ and the results for the velocity- 
dependence of noise intensities [lj} to give a tentative 
explanation of this phenomenon. 

The additive noise intensity considerably increases for 
bigger \u\ due to the relationship a ~ e' u ' (incoher- 
ent noisy background responsible for the random walk 
behavior of acceleration is much intensified) while the 
multiplicative noise intensity remains at approximately 
the same level (very intense vortical structures responsi- 
ble for the random multiplicative process are relatively 
frozen in time and characterized by a saturation level 
of \u\ for a given Reynolds number and vorticity). This 
means an increase of large scale effects produced by the 
interaction between small and large scales (the effect of 



7 



nonlocal interactions) . Since the large scales are the only 
unifying agent between the noises the cross correlation 
of the noises becomes more pronounced for bigger \u\. 
Whereas a direct effect of the increase of additive noise 
intensity tends to symmetrize the distribution (see Fig.[SJ) 
due to a higher degree of chaoticity (higher statistical 
isotropy), bigger |u| may imply a variation of the cross 
correlation, to which the skewness is highly sensitive, so 
that the overall effect is a small but appreciable increase 
of the mean of the conditional distribution. To be more 
precise, the net effect depends on competition between A 
and a as the velocity |u| increases. 

Also, it is worthwhile to note that the effect of dis- 
crete Kolmogorov turbulent cascade may be of impor- 
tance here since it is characterized by a relationship be- 
tween high-amplitude harmonics of the basic ratio (non- 
local interactions). 

The following remark is in order. 

Note that in Sec. Ill Bl we used the experimental data 
on unconditional P(a) instead of that for the conditional 
P(a\u) with the aim to estimate result of the possible 
dependence of A on u. Such an approach is fully justified 
since these curves are of a similar character and we were 
interested only in the general effect of the variation of 
A = X(u) on the shape of P(a\u), for which case partic- 
ular values of the fitted parameters are not important. 
This dependence was found to be not capable to explain 
the characteristic variation of the shape of P(a\u) with 
variation of u shown in Fig. 6a of Ref. [l^. I n other 
words, it does not make considerable contribution to the 
variance of P(a\u) consistent with the experimental data. 
However, this evidently does not exclude possible pres- 
ence of a specific, e.g., polynomial or exponential, depen- 
dence of A on u which would allow one to qualitatively 
explain the observed small variation of the conditional 
mean acceleration {a\u} with variation of u. 

Note that we have evaluated the conditional mean 
acceleration normalized to unit conditional acceleration 
variance (a 2 !^} 1 / 2 rather than to unit unconditional ac- 
celeration variance (a 2 ) 1 / 2 . This makes a difference since 
the conditional acceleration variance (a 2 ^) 1 / 2 increases 
with an increase of u 0, contrary to the case of 
unconditional acceleration variance that is constant. 

In the RIN approach, the evaluation of uncondi- 
tional acceleration variance requires a calculation of the 
stochastic expectation of P(a\u) over random u. Whereas 
we can judiciously choose, in the first approximation, the 
distribution of u to be Gaussian with zero mean and 
ct = e u / u ° jToL Il3l ] we should also take into account for 
the emerging effect of A, for which case we have no clue 
to make a judicious choice of some function X(u). An 
accounting for A(u) may produce a considerable contri- 
bution to P(a) due to the integral effect. For this reason, 
in the present paper we do not evaluate P(a) unless we 
identify the functional form of X(u). This issue is of much 
interest and can be studied elsewhere. 

Nevertheless, one can proceed here with a model- 
independent study by normalizing the conditional mean 



acceleration to unit conditional acceleration variance at 
u = which is constant (does not depend on u), i.e., by 
evaluating {a\u) / (a^O) 1 / 2 . The result is plotted in Fig. [5] 
(boxes) and shows a much steeper decrease as compared 
with that of {a\u) / (a 2 ^) 1 / 2 (triangles), but yet no in- 
crease. 

In summary, the predicted conditional mean acceler- 
ation (a|u)/(a 2 |0) 1 ^ 2 is negligibly small and decreases 
down to zero for larger |u| under the assumption that 
a _ e l u |/3 an( j -(. ne othgj. parameters to be constant. 
Hence, the additive noise tends to symmetrize acceler- 
ation distribution. 

Therefore, we are left with the only possibility to ex- 
plain the observed appreciable increase of (a\u) / (a 2 ) 1 / 2 
with the increase of \u\: to assign some dependence of the 
parameter A, which measures correlation between the ad- 
ditive noise and the multiplicative one, on velocity fluc- 
tuations u. It should be noted that this is in agreement 
with the LDN model, in which the additive noise depends 
on small-scale velocity fluctuations. 

It should be stressed however that as it has been men- 
tioned above the observed nonzero mean acceleration is 
associated to the flow anisotropy which may be not re- 
lated to the effect described by the parameter A. 



IV. DISCUSSION AND CONCLUSIONS 

(i) We have shown that the cross correlation parameter 
A of the LDN type model Q for a particle acceleration 
could be used to explain a skewness of the acceleration 
probability density function, and estimated its value, A ~ 
—0.005, by using a fit to the recent experimental statistics 
data on the transverse component of acceleration. 

(ii) The mean acceleration is found to be very close 
to zero, (a|u)/(a 2 |0) 1/2 < 0.0024, when the predicted 
a A P(a) is fitted to the experimental data. The mean ac- 
celeration should vanish for homogeneous isotropic tur- 
bulence. The observed mean acceleration can be at- 
tributed to small anisotropy (imperfection) of the studied 
R\ = 690 flow. However, the DNS, for which isotropy 
is well satisfied, indicates slight departures from zero. 
Whether this is of some importance for developed turbu- 
lence should be clarified. In any case, the observed no- 
ticeable nonzero mean of Lagrangian acceleration, which 
is usually associated to extremely small scales of the flow, 
in the high-Rcynolds-numbcr flow that is anisotropic at 
large scale, deserves a separate study. 

(iii) Using the RIN approach which extends the LDN 
type model by assuming certain relationship of noise in- 
tensities and in general other model parameters to ve- 
locity fluctuations u we have studied acceleration statis- 
tics conditional on velocity fluctuations. We found that 
the assumption A = A(u) could not be responsible for 
the experimentally observed characteristic variation of 
the shape of conditional distribution P(a\u) with varia- 
tion of u. Taken together with the result of our previous 
work |Toj this implies that only the additive noise inten- 
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sity a reveals an essential dependence on u for which we 
used the exponential function, a — e u ' u °, relevant from 
both the phenomenological and experimental points of 
view. The additive noise tends to symmetrize accelera- 
tion distribution for larger \u\: the simulated conditional 
mean acceleration is very small and decreases down to 
zero as shown in Fig. [S] This is not in agreement with 
the experimental data even qualitatively. However, the 
observed mean acceleration could be attributed mainly to 
anisotropy (imperfection) of the studied R\ — 690 flow. 

(iv) We have found that despite the previously ob- 
tained result that a sole dependence of the additive 
noise intensity parameter a on it is in a good qualitative 
agreement with the experimental data on the conditional 
acceleration variance, (a 2 \u) j (a 2 ) , it appears to be not 
capable to capture qualitatively the experimental data on 
the conditional mean acceleration, (a|w)/(a 2 ) 1 / 2 , which 
increases for bigger \u\. While determination of the con- 
ditional mean acceleration for higher isotropic turbulent 
flows can be left for future experiments and numerical 
simulations, we speculate on the relationship between A 
and skewness. This suggests that, in addition to the de- 
pendence a = a(u), the cross correlation parameter A, 
which directly controls the skewness, could depend on u 
in some way. 

(v) As the result, to meet the available experimen- 
tal data and the DNS we are led to consider, in 
a self-consistent RIN approach, the conditional prob- 
ability density function © in the form P{a\u) — 
P(a\a(u), A(w)), with a — e u l u °, velocity fluctuations u 
to be Gaussian distributed with zero mean, some func- 
tion A(it) to be determined, and free parameters D, B, 
and v$ (not depending on u, in the first approximation), 
to be used for a fitting to the experimental P(a\u). This 
study is of interest and can be made elsewhere. 

(vi) The marginal distribution P m (a) is evaluated by 
integrating out u, 

/oc 
duP(a\a(u),\(u))g(u), (10) 
-oo 

where in the first approximation the distribution g{u) can 
be taken Gaussian, and then should be compared with 
the experimental unconditional distribution. This re- 
quires prior determination of the dependence \{u) which 
can be made elsewhere. 

We conclude by a few remarks. 

(a) As it has been mentioned above, the distribution 
© can be used "as is" for a fitting to the experimen- 
tal unconditional P(a), as shown in Fig. 2] It should 



be stressed however that in the proper RIN approach 
it is the marginal distribution l|10(l obtained from © 
by averaging over random noise intensities with some 
judiciously chosen distributions assigned to them, that 
should be fitted to the experimental unconditional P(a). 
The simplest choice is to assume that only a is a ran- 
dom parameter, inverse of which follows chi-square or 
log-normal distribution in the spirit of simple RIN mod- 
els 01 • Particularly, for a = e u / u ° the choice of the log- 
normal distribution is equivalent to that u is normally 
distributed with zero mean It should be empha- 

sized here that only absolute value of u contributes to the 
marginal distribution. In this case the distribution J3J) is 
treated as a conditional probability distribution function 
P{a\a(u)) (see Fig. 7 in Ref. (Ty|), which could be in 
principle fitted to the experimental conditional distribu- 
tion P{a\u). Although we have found a good qualitative 
agreement of the model with the experimental data on 
conditional statistics, an illustrative qualit y o f the rep- 
resentation of experimental P(a\u) in Ref. [ljj does not 
allow us to make a reliable numerical fit of the proposed 
P(a\u) — P(a\a(u), X(u)). For good fit results, high ac- 
curacy experimental data on P(a\u) and on the contri- 
bution to fourth order conditional moment, a A P{a\u), for 
| it | ranging from zero to three rms velocity with the step 
0.5, would be required. 

(b) This would give a possibility to identify the de- 
pendence A(it), for which one can try a polynomial or 
exponential function. 

(c) Despite similarity one observes some differ- 
ence between the experimental distributions P(a) and 
P(a\u) 19] which could be described in terms of the RIN 
approach to LDN type model © along the line of rea- 
soning given in the present paper. Namely, with P(a\u) 
taken to be the distribution P(a\a(u), A(u)) and P(a) to 
be derived from it by integrating out velocity fluctua- 
tions, Eq. [flfjjl . 

(d) Reynolds-number dependence of the acceleration 
statistics is not considered in the present paper and can 
be studied elsewhere. Also, it is of interest to study vari- 
ation of the value of flatness factor of acceleration distri- 
bution measuring Lagrangian intermittency as a function 
of finite time increment used for low-pass filtering. This 
dependence exhibits a fine structure of the viscous dis- 
sipation range of time-scales. The flatness factor of the 
R\ = 690 flow varies by about 15% when the filter width 
is changed from 0.23^ to 0.31r^, where t v is Kolmogorov 
time [lj|. 
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